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ABSTRACT 

We study the globular cluster (GC) system of the dust-lane elliptical galaxy NGC 6702, 
using B, V and /-band photometric observations carried out at the Keck telescope. 
This galaxy has a spectroscopic age of « 2 Gyrs suggesting recent star-formation. 
We find strong evidence for a bimodal GC colour distribution, with the blue peak 
having a colour similar to that of the Galactic halo GCs. Assuming that the blue 
GCs are indeed old and metal-poor, we estimate an age of 2-5 Gyrs and supersolar 
metallicity for the red GC subpopulation. Despite the large uncertainties, this is in 
reasonable agreement with the spectroscopic galaxy age. Additionally, we estimate a 
specific frequency of Sn = 2.3 ± 1.1 for NGC 6702. We predict that passive evolution 
of NGC 6702 will further increase its specific frequency to Sn ~ 2.7 within 10 Gyrs, 
in closer agreement to that of typical present-day ellipticals. We also discuss evidence 
that the merger/accretion event that took place a few Gyrs ago involved a high gas 
fraction. 

Key words: globular clusters: general - galaxies: individual: NGC 6702 - galaxies: 
star clusters. 



; 1 INTRODUCTION 



The 'merger' hypothesis, postulating that elliptical galax- 
ies form by spiral galaxy mergers, was first proposed by 
Toomre & Toomre (1972). This scenario has gained new 
impetus over the last few years by both numerical simula- 
tions (e.g. Barnes 1992; Hernquist 1992; Mihos & Hernquist 
1996) and observational studies, suggesting that the prop- 
erties of merger remnants are consistent with those of el- 
lipticals (e.g. Toomre 1977; Schweizer et al. 1990; Casoli et 
al. 1991; Schweizer & Seitzer 1992; Hibbard & van Gorkom 
1996; Georgakakis, Forbes & Norris 2000). 

However, one of the long standing drawbacks of the 
'merger' hypothesis is the excess globular clusters (GCs) per 
unit starlight around ellipticals compared to spirals (van den 
Bergh 1984). To overcome this problem Schweizer (1987) 
and Ashman & Zepf (1992) suggested that the starburst ac- 
tivity usually experienced during gas-rich mergers will form 
new GCs and thus increase their number compared to their 
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progenitors. Indeed, there is accumulating observational evi- 
dence supporting this scenario in on-going mergers (e.g. Lutz 
1991; Whitmore et al. 1993; Whitmore & Schweizer 1995; 
Schweizer 1996; Miller et al. 1997; Zepf et al. 1999; Forbes & 
Hau 2000). However, the situation in ellitpical galaxies (the 
eventual product of mergers) is less clear. Although many 
ellipticals have been observed to have bimodal colour distri- 
butions, when examined in detail their GC properties do not 
match those from the merger picture (Forbes et al. 1997). 
Outstanding questions include whether the newly formed 
GCs survive and whether the final GC system resembles that 
of typical (old) ellipticals. Probing the GC properties in in- 
termediate age merger-remnants (1-5 Gyrs) may be of key 
importance for addressing these issues. Indeed, such young 
ellipticals may provide the link between on-going mergers 
and relaxed ellipticals. 

For example, a sample of merger-remnants spanning a 
range of ages will allow study of the chronological devel- 
opment of the system's specific frequency (measuring the 
number of globular clusters per unit starlight) in compari- 
son with that of evolved ellipticals. Young ellipticals are also 
expected to have two distinct GC populations with differ- 
ent colour and magnitude distributions. By exploiting these 
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differences one can attempt to age-date the newly formed 
GCs, which in turn can be compared to the time since the 
gaseous merger. Indeed, studies of the dynamically young 
proto-ellipticals NGC 1700 and NGC 3610 have shown the 
presence of bimodality in their colour distributions (Whit- 
more et al. 1997; Brown et al. 2000). Also the age of the 
new GCs, although model dependent, is found to be similar 
to estimates for the galaxy age based on spectroscopic or 
dynamical arguments. 

In this paper we study the globular cluster system of 
the elliptical galaxy NGC 6702. Lauer (1985) found a dust 
lane in the centre of NGC 6702 indicating a recent gaseous 
merger event. This galaxy is also relatively luminous at far- 
infrared wavelengths compared to other elliptical galaxies, 
with Lfir = 1.3 x 10 10 L (Wiklind et al. 1995). Terlevich 
& Forbes (2000) have estimated a spectroscopic age and 
metallicity of about 2Gyrs and [Fe/H] = +0.5 respectively 
for NGC 6702 suggesting that it is a young elliptical. The 
main observational properties of NGC 6702 are given in Ta- 
ble 0. 

Additionally, according to the NED database the near- 
est galaxy projected on the sky that is either brighter 
or no more than 2 magnitudes fainter, than NGC 6702 is 
NGC 6703 at a projected distance of 10 arcmin. However, 
as NGC 6702 has a recession velocity of 4727 km s -1 and 
NGC 6703 has V — 2365 km s -1 they are actually separated 
by some 30Mpc assuming no peculiar motions. The near- 
est galaxy with a similar recession velocity is NGC 6711 
(V = 4671 kms -1 ) with a separation of ~ 120 arcmin or over 
2 Mpc. Other galaxies listed by NED are more than 2 magni- 
tudes fainter than NGC 6702. Thus NGC 6702 is an example 
of an isolated elliptical galaxy with no comparable galaxies 
within 2 Mpc. Such galaxies are fairly rare and hence their 
evolutionary history is of some interest. In particular, ellip- 
ticals in low density environments are interesting within the 
hierarchical merger framework as they are expected to have 
assembled much more recently than their cluster counter- 
parts (Kauffmann & Chariot 1998). Any merger would be 
expected to have occurred in the more recent past and in- 
volved progenitors that have largely processed gas into stars. 

In section |i| the observations, data reduction and the 
isophotal properties of NGC 6702 are presented. Section 
^| describes the candidate GC selection, while section ^ 
presents our results which are discussed in section [| Finally, 
section ^ summarises our conclusions. Throughout this pa- 
per we adopt Ho = 75 kms -1 Mpc -1 . 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Description of the data 

Broad-band imaging of NGC 6702 in the B, V and I filters 
was carried out at the Keck-II telescope on 1999 August 17, 
using the the Low Resolution Imaging Spectrometer (LRIS; 
Oke et al. 1995). The LRIS instrument, equipped with a 
TEK 2048 x 2048 CCD, is mounted on the Cassegrain focus 
providing a 0.215 arcsec pixel -1 imaging scale and a 6' x 8' 
field-of-view. The total exposure time in the B, V and / 
filters was 3600, 1800, and 1200 sec respectively, segregated 
into three separate integrations of 1200, 600 and 400 sec re- 
spectively. The seeing at the time of the observations varied 
between 0.8 to 1 arcsec. 




Figure 1. V-band image of NGC 6702 obtained at the Keck tele- 
scope. 

The data were reduced following standard procedures, 
using IRAF tasks. The reduced images are found to be flat 
to better than ~ 2%. The final V-band image of NGC 6702 
is shown in Figure |l[ Photometric calibration was performed 
using standard stars from Landolt (1992). The photometric 
accuracy, estimated using these standard stars, is ±0.02 mag 
in all three bands. We adopt a Galactic extinction in the 
NGC 6702 direction of E(B—V) = 0.06, taken from Burstein 
& Heiles (1984). Assuming a standard reddening curve (Sav- 
age & Mathis 1979), we correct each magnitude for Galac- 
tic extinction using the values Ab = 0.25, Ay = 0.19 and 
Ai = 0.09 mag. 

2.2 Galaxy surface brightness profile 

The _B-band surface brightness profile, /is, of NGC 6702 is 
plotted in Figure [| This has been determined by fitting el- 
lipses to the galaxy intensity profile using the ISOPHOTE 
task within STSDAS. During the ellipse fitting process the 
centre of the galaxy was kept fixed and a 3 a clipping al- 
gorithm was employed. Also, bright objects as well as the 
central saturated parts of the galaxy, were masked out. The 
position angle (PA) and ellipticity (e) were fitted at each 
radius, until S/N constraints terminated the fitting process 
at low surface brightness. For NGC 6702 this occurs at an 
equivalent radius r eq « 50" (r eq = a y/l — e, where a is the 
semi-major axis). Beyond that radius the position angle and 
ellipticity were kept constant. The sky level was determined 
using the method described by Goudfrooij et al. (1994). 

Although the isophotes of elliptical galaxies are well ap- 
proximated by ellipses, deviations at the few per cent level 
are common. Such deviations are often quantified by the 4th 
sine and cosine amplitudes (S4 and C4 respectively) of the 
Fourier transform of the galaxy intensity at a given radius. 
The standard interpretation of the residual C4 terms is that 
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Name RA (J2000) DEC (J2000) 1 b type heliocentric distance B°, V% 

(RC3) velocity (kms -1 ) modulus (mag) (mag) 
NGC6702 18 h 46 m 57.6 s +45° 42' 20" 75.03° 19.79° E3 4727 34.00 13.04 12.15 



Table 1. Observational properties of NGC 6702 obtained from NED. The distance modulus is estimated assuming Hq = 75 km s 1 Mpc 



they indicate either the presence of a weak underlying disk 
if positive or boxiness if negative. However, it has been em- 
phasised that both C4 and S4 terms should be considered 
when estimating the isophotal shape of galaxies (Franx et 
al. 1989; Peletier et al. 1990; Goodfrooij et al. 1994). For ex- 
ample, the presence of a disk, bar or ring with a significant 
position angle twist relative to the elliptical galaxy semi- 
major axis is expected to give both positive and negative 
C4 and/or S4 terms, depending on its orientation. Conse- 
quently, in the case of such projection effects we expect S4 
terms to be significant. 

Figure ^ plots the ellipticity, position angle as well as 
C4 and S4 terms as a function of radius for NGC 6702 in 
the B-band. Similar trends are apparent for the V and I- 
bands. The C4 term is negative at small galactocentric dis- 
tances, indicating boxiness, i.e. excess light at 45° relative 
to the semi-major axis. This is partly due to the presence 
of the dust lane in the core of NGC 6702 (Lauer 1985) lying 
perpendicular to the semi-major axis and extending out to 
~ 20". The S4 term is also non-zero within the same range of 
r eq , albeit with a smaller amplitude. Masking out the dust 
lane reduces the amplitude of both the C4 and S4 terms. 
These terms are also non-zero (S4 positive and C4 nega- 



tive) at 



40 indicating boxiness. At larger distances 



(r eq > 60") both the C4 and S4 deviate significantly from 
the elliptical symmetry, but this is likely to be an artifact 
arising from the low S/N at these radii. 




r eq (arcsec) 

Figure 2. B-band surface brightness profile of NGC 6702. The 
dashed lines are the la uncertainty envelopes. We define the 
equivalent radius r eq = a \J1 — e, where a and e are the semi- 
major axis and ellipticity respectively. 



2.3 Source extraction and completeness limits 

To detect sources superimposed on NGC 6702 we subtract 
the galaxy light profile using a median filter technique. 
Firstly, the reduced image for each filter is smoothed us- 
ing a median rectangular filter with a size of 40 pixels. The 
smoothed image is then subtracted from the original frame 
and the resulting image is subsequently used for source ex- 
traction and photometry. 

The sources in the reduced frames were extracted us- 
ing the SExtractor package (version 2.1.0; Bertin & Arnouts 
1996). The main input parameters are the detection thresh- 
old, given as a multiple of the sky variance (asky) and the 
minimum area in pixels for an object to be extracted. Af- 
ter considerable experimentation, we adopted a threshold 
of 2.5x<rsfc y and a minimum area of 10-15 pixels (~0.5- 
0.7 arcsec 2 ). This choice of values minimises the number 
of spurious detections, while ensuring that faint objects 
are successfully extracted. A total of 1455, 1456 and 1494 
sources were extracted form the B, V and I frames respec- 
tively. The photometric catalogues generated for each wave- 
band were then matched, resulting in a total of 939 common 
detections in all three filters, B, V and /. This source list is 
referred to as the matched catalogue. Using this catalogue 
we performed photometry using the PHOT task. After a 



curve-of-growth analysis on several sources an aperture of 
8 pixels radius was employed with the background estimated 
from an annulus of 15 to 20 pixels around each source. For 
unresolved sources the 8 pixels radius aperture magnitude 
is close to total. 

The completeness limit in a given filter is estimated by 
adding artificial point sources of different magnitudes to the 
original image, using MKOBJECTS task. Then we attempt 
to recover the artificial objects using SExtractor package, 
with the same parameter settings as with those adopted for 
the real data. A total of 1500 point sources were added to 
a given filter frame. To avoid overcrowding the frame, ar- 
tificial sources are added in groups of 150 and the process 
is repeated 10 times. The results are presented in Figure ^j, 
plotting the fraction of recovered objects as a function of 
magnitude in the B, V and I filters respectively. The 80% 
completeness limits in the B, V and I filters are found to 
be ~ 26.1, ~ 25.3 and « 23.9 mag respectively. The mean 
photometric error at these magnitude limits is « 0.1 mag for 
all three filters. We have also investigated variations of the 
completeness limits due to the increased background close to 
the galaxy centre. A method similar to that described above 
was adopted, with artificial point sources being added within 
successive elliptical rings centred on NGC 6702 and having 
the same mean ellipticity and position angle as the galaxy 
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Figure 3. Radial profile in the B-band of the ellipticity (e), po- 
sition angle (PA), 4th cosine term (C4) and 4th sine term (S4), 
estimated by fitting ellipses to the galaxy light profile. At radii 
larger than 50" both the ellipticity and the position angle are 
kept constant during the ellipse fitting routine, due to S/N con- 
straints. 



Figure 4. Fraction of successfully recovered point sources as a 
function of magnitude in the B (dashed line), V (continuous line) 
and I (dotted line) filters. The vertical lines represent the peak 
magnitude of the Milky Way GC luminosity function in different 
wavebands (Ashman & Conti 1995; Mj iPea t = —6.50, Mv, pe .ak = 
-7.33, M I:Peak = -7.95), shifted to the distance of NGC6702. 



(e = 0.2, PA=145°). The effect of galaxy background was 
found to be small except for r £q < 25". However, this is close 
to the saturated region of the CCD images which is masked 
out in our analysis. Therefore, the lower sensitivity close to 
the galaxy center will not affect the results from this study. 



3 SELECTION OF GLOBULAR CLUSTER 
CANDIDATES 

Globular cluster candidates are selected from the matched 
catalogue on the basis of their expected range in colour and 
magnitude. In particular, the full range of possible glob- 
ular cluster metallicities is —2.5 < [Fe/H] < +1.0. Us- 
ing the colour-met allicity relation of Couture et al. (1990) 
and Kissler-Patig et al. (1998), the above mentioned metal- 
licity interval corresponds to 1.2 < B — / < 2.5 and 
0.6 < V - I < 1.7. Figures | and | plot the B — I and 
V — I colours against B — V respectively, for the objects 
in the matched catalogue. The horizontal lines in these fig- 
ures represent the above mentioned B — I and V — / colour 
cutoffs. Also shown are the regions of the parameter space 
occupied by the Galactic GCs (Harris 1996; corrected for 
Galactic reddening). The accepted range of B — V colour 
for globular cluster selection is chosen to span the B — V 
colour range of the Galactic GCs. Additionally, using the 
population synthesis models of Worthey (1994) we find that 
the B — I, V — I and B — V colours of a stellar population 
of a given age get redder with increasing metallicity in the 
direction indicated by the arrows in Figures |H] and q. 

Globular cluster candidates are selected on the basis of 



their photometric colours if they lie within the regions indi- 
cated in Figures ^ and |(| adopting a 3<r uncertainty on each 
colour given by the PHOT software. A total of 374 candidate 
GCs are selected. These objects are then visually inspected 
to discard resolved sources, likely to be galaxies, and sources 
lying near hot pixels. The final list consists of 205 candidate 
GCs. Although visual inspection is not an objective method 
for identifying point sources, we believe that to the first ap- 
proximation we successfully eliminate background galaxies 
from the GC catalogue. This is demonstrated in Figure ^ 
plotting an estimator of the spatial extent of a source as a 
function of /-band magnitude. The size of a source is esti- 
mated by the difference between the /-band magnitude mea- 
sured within 3 and 8 arcsec radius apertures. It is clear from 
Figure (7] that visually classified resolved sources are sepa- 
rated from point-like objects, particularly at bright magni- 
tudes. At fainter magnitudes (/ > 23 mag), however, the 
visual classification scheme becomes more uncertain. Also 
note that point sources brighter than / ~ 20 mag are likely 
to have saturated cores. 

Additionally, we apply a bright magnitude cutoff set 
to ~ 2 mag brighter than the most luminous Milky Way 
GC (fl-Cen; M v = -10.29 mag), shifted to the distance 
of NGC 6702. This is because NGC 6702 has a young stellar 
component, with an estimated spectroscopic age of « 2 Gyrs 
(Terlevich & Forbes 2000). Therefore, NGC 6702 may pos- 
sess a number of young metal-rich GCs, expected to be sig- 
nificantly brighter than the old metal-poor GCs. For ex- 
ample, the Worthey (1994) models predict a 2.0 Gyrs old 
stellar population to be brighter by as much as 1.5 mag in 
the /-band compared to an old metal-poor stellar popula- 
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Figure 5. B — I , B — V colour-colour diagram for the sources 
detected in all three niters, B, V and / (dots). The Milky Way 
globular clusters are shown as filled circles. The lines delimit the 
region of the parameter space that globular clusters are expected 
to occupy (see text for details). The arrow indicates the direction 
the B—V and B — I colours redden for a single burst stellar popu- 
lation model of constant age and increasing metallicity (Worthcy 
1994). 



Figure 6. V — I, B — V colour-colour diagram for the sources 
detected in all three filters, B, V and / (dots). The Milky Way 
globular clusters are shown as filled circles. The lines delimit the 
region of the parameter space that globular clusters are expected 
to occupy (see text for details). The arrow indicates the direction 
the B — V and V — I colours redden for a single burst stellar popu- 
lation model of constant age and increasing metallicity (Worthey 
1994). 



tion (15Gyrs, [Fe/H] = —1.5) assumed to be representative 
of the Milky Way GCs. The bright magnitude cutoffs are 
thus set to B — 22.5, V = 21.5 and / = 20.5 mag. Sources 
brighter than this are almost certainly stars. In the rest of 
the paper, unless otherwise stated, we consider GCs brighter 
than the 80% completeness limit for point sources. There- 
fore, the magnitude range of the candidate GC sample is 
22.5 < B < 26.1, 21.5 < V < 25.3 and 20.5 < I < 23.9 mag. 
The colour and magnitude selected and visually inspected 
sample comprises a total of 151 GC candidates. It should 
be noted that including GCs brighter than the 50% com- 
pleteness limit for point sources does not change any of our 
conclusions. 

Figure ^ shows the colour-magnitude diagram of all the 
colour selected and visually inspected GC candidates before 
applying any magnitude cuts. Objects have a wide range of 
colours but there is some evidence for two concentrations 
at B - I « 1.5 and B - I « 2.3 mag. This will be further 
discussed in section ^. 

Although we attempt to minimise the fraction of fore- 
ground stars and background galaxies in the final GC sam- 
ple, it is likely that a number of contaminating sources are 
also selected. 

The number counts and colour distribution of fore- 
ground stars are estimated using the model developed by 



Bahcall & Soneira (1980)^. In this model the Galaxy con- 
sists of both a disk and a spheroid component from which 
the distribution of stars in the B and IZ-bands as well as the 
B — V colour distribution are calculated. To predict both 
the star counts in the /-band and the B — I and V — I 
colour distributions we use the empirical tight relation be- 
tween the B — V and V — I colours for stars (Bessel 1990). 
We predict about 120 stars within the field-of-view in the 
magnitude range 22.5 < B < 26.0, 21.5 < V < 25.5 and 
20.5 < I < 24.0 mag. However, the colour selection (assum- 
ing a 3ct uncertainty of 0.2 mags) eliminates the majority of 
stars from the final catalogue. Indeed, the model of Bahcall 
& Soneira (1980), after both magnitude and colour selection, 
predicts a total of about w 20 stars in the final GC sample 
(total of 151 sources). Therefore, the star contamination is 
about 13%. 

To quantify the contamination of the final GC sample 
by background galaxies we should first estimate the com- 
pleteness limits for galaxies at different wavebands. Since 
galaxies are resolved sources, their completeness limit is 
expected to be at brighter magnitudes to that of point 
sources (e.g. GCs in NGC 6702). Figure ^| plots the num- 
ber counts for all the sources detected in our B-band frame 
(total of 1455). Also shown is the star count prediction from 
the model of Bahcall & Soneira (1980). The difference be- 



3 The code for the star cou nt model can be obtained from 
bttp:/ /ascl.net/bsgmodel.html 
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Figure 7. Size parameter against /-band magnitude. The size 
parameter, estimating the extent of a source, is defined as the I- 
band magnitude difference within 3 and 8 arcsec radius apertures, 
ft is clear that visually classified resolved sources (filled circles) 
are reasonably well separated from point-sources (open circles), 
although there is increasing scatter at fainter magnitudes. 
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Figure 8. Colour— magnitude diagram for the 205 candidate GCs, 
after colour selection and visual inspection. No magnitude cutoff 
has been applied. The continuous line corresponds to the com- 
pleteness limit / = 23.9 mag. The dashed lines are the B-band 
magnitude cutoffs for the selection of the final candidate GC sam- 
ple. 



tween the total observed source counts and the predicted 
star counts in a given magnitude bin estimates, to the first 
approximation, the number of background galaxies detected 
in our field (here we ignore the contribution from GCs to 
the counts). Also shown in Figure |^ are the B-band galaxy 
counts from Metcalfe et al. (1991, 1995). It is clear that the 
predicted galaxy counts in our field are in good agreement 
with those from Metcalfe et al. (1991, 1995) to the magni- 
tude B m 25 and then decline due to incompleteness. We 
estimate our B-band galaxy catalogue to be about S3 50% 
and S3 30% complete at S3 25.5 and S3 26.0 mag respectively. 
Figures and are similar to Figure |^ for the V and 
/-bands respectively. From these Figures, we estimate our 

V and /-band galaxy catalogues to be S3 50% complete at 

V S3 24.5 and / S3 23.5 mag respectively. 

Ignoring GCs when calculating the number of galaxies 
within the field-of-view at B, V and /-bands, overestimates 
the galaxy numbers at faint magnitudes. Consequently, the 
above mentioned completeness limits are conservative esti- 
mates. 

To quantify the contamination of the GC sample due to 
background galaxies we use the multiwavelength (UBVI) 
faint galaxy catalogue of Lilly, Cowie & Gardner (1991) 
complete to Iab = 24.5. The AB system is defined as 
AB = 48.6 — 2.5 log /„, where /„ is the flux cgs units. The 
transformations from the AB system are: B S3 Bab + 0.17, 

V S3 Vab and / S3 l AB - 0.48 (Lilly, Cowie & Gard- 
ner 1991). Using this galaxy sample we predict a total of 
S3 530 galaxies with 22.5 < B < 26, 21.5 < V < 25.5 
and 20.5 < I < 24 mag after converting to our magnitude 
system and scaling to the area of the LRIS field-of-view 
(S3 31 arcmin 2 ). Taking into account the galaxy incomplete- 
ness of our sample, estimated from Figures l|, we predict 
a total of ss 320 galaxies within the above mentioned mag- 
nitude limits. Colour selection, using the aperture colours 
of Lilly, Cowie & Gardner (1991) and assuming a 3a un- 
certainty of 0.2 mags, reduces the number of contaminating 
background galaxies in our candidate GC sample to w 160. 

However, as demonstrated in Figure |^, visual inspection 
is efficient in removing extended sources from the sample 
especially at bright magnitudes. Indeed, a total of w 120 
resolved sources with 22.5 < B < 26, 21.5 < V < 25.5 and 
20.5 < / < 24.0 mag have been identified, leaving about 
ss 40 unidentified galaxies within the final GC sample. This 
corresponds to a contamination of the GC candidate cata- 
logue (total of 151) by galaxies of < 25%. We stress that this 
fraction is an upper limit. In particular, we simulated galaxy 
bulges following r 1//4 -law profile, with absolute luminosity 
and effective radius of Ms=-21.0mag and r e ff = 4kpc re- 
spectively, typical of L* ellipticals. We found that only half 
of them are detected to the limit B — 25.5 mag and this 
fraction drops to 10% at B — 26.0 mag. Moreover, disk- 
dominated galaxies will be missed by the detection algo- 
rithm at even brighter magnitudes. These surface bright- 
ness effects are likely to be significantly less severe in the 
deep surveys used to determine the galaxy counts at faint 
magnitudes. Indeed these surveys have effective exposure 
times as much as 10 times longer than those used here and 
are therefore, specifically designed to detect sources at faint 
magnitudes. Since bulges and compact galaxies account for 
less than 20% of the galaxy population at faint magnitude 
limits (Elson & Santiago 1996), we anticipate the galaxy 
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Figure 9. Differential number counts against B-band magnitude 
for all the sources detected in the NGC 6702 B-band frame (his- 
togram). The dashed line is the star counts predicted by the model 
of Bahcall & Soneira (1980). The open circles are the estimated 
number of galaxies on the NGC 6702 B-band frame (see text). 
Filled circles are the B-band galaxy counts of Metcalfe et al. 
(1991, 1995). It is clear that at bright magnitudes the estimated 
number of galaxies detected on the NGC 6702 B-band frame are 
in fair agreement with those predicted by Metcalfe et al. (1991, 
1995). At faint magnitudes (B > 25.0 mag) incompleteness is af- 
fecting the present galaxy sample. 



Figure 10. Differential number counts against V-band magni- 
tude for all the sources detected in the NGC 6702 y-band frame 
(histogram). The dashed line is the star counts predicted by the 
model of Bahcall & Soneira (1980). The open circles are the 
estimated number of galaxies on the NGC 6702 V-band frame 
(see text). The continuous line is the best fit to the faint V- 
band galaxy counts estimated by Woods & Fahlman (1997). At 
bright magnitudes the estimated number of galaxies detected on 
the NGC 6702 V-band frame are in fair agreement with those 
predicted by Woods & Fahlman (1997). At faint magnitudes 
(V > 23.0 mag) incompleteness is affecting the galaxy counts from 
our study. 



contamination of our final candidate GC sample to be much 
less than ~ 25%. 



4 RESULTS 

4.1 Colour distributions 

The B — I colour distribution for the final GC sample is plot- 
ted in Figure There is evidence for bimodality with the 
peaks of the distribution at B — / « 1.50 and B — I ~ 2.15. 
Indeed the KMM statistical test (Ashman et al. 1994) re- 
jects the single Gaussian model at a confidence level better 
than 99.99%, giving peaks at B — I = 1.52 and 2.16 respec- 
tively. It should be noted that this result is obtained after 
clipping objects with B — I > 2.6 mag, likely to be back- 
ground galaxies (3 out of 151; see Figure |l3|) . However, this 
is not expected to artificially increase the likelihood of the 
two Gaussian fit. In particular, Ashman et al. (1994) argue 
that clipping of the data before applying the KMM test is 
sometimes desirable to remove background contamination. 
Also shown in Figure ^ is the reddening-corrected colour 
distribution of Galactic GCs that comprise old metal-poor 
GCs. The blue peak of NGC 6702 GCs lies close to the Milky 
Way GC distribution. 

Figure [l3| compares the candidate GC B — I colour dis- 
tribution with the colours of contaminating stars and galax- 



ies (after being selected on the basis of their magnitudes 
and colours) predicted in the previous section. It is clear 
that foreground stars do not significantly modify the ob- 
served candidate GC colour distribution. Background galax- 
ies, however, have a larger effect on the observed colour 
distribution, having colours similar to those of the red GC 
peak. Nevertheless, the fraction of contaminating galaxies 
within the sample is expected to be sufficiently small that it 
will not significantly modify the position of the peaks of the 
GC colour distribution. We further investigate this by sub- 
tracting the predicted galaxy colour distribution from that 
of the GC candidates in Figure [H| The KMM statistical 
test is then applied to the resulting distribution. A single 
Gaussian model is rejected at a confidence level better than 
99.99%, while the peaks of the bimodal distribution are at 
B — I — 1.56 and 2.27 mag respectively. Although both peaks 
are redder than those recovered using the full candidate GC 
colour distribution, the difference is sufficiently small that 
does not alter any of our conclusions. 

Figure |l4| plots the V — I colour histogram of both the 
NGC 6702 GC candidates and the Galactic GCs. There is no 
strong evidence for bimodality, with the KMM test finding 
no improvement of the two-group fit over a single Gaussian. 
This seems surprising since the B — I colour distribution 
shows clear evidence for bimodality. It is likely that the ab- 
sence of bimodality in Figure O is due to the smaller wave- 
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Figure 11. Differential number counts against /-band magnitude 
for all the sources detected in the NGC 6702 /-band frame (his- 
togram). The dashed line is the star counts predicted by the model 
of Bahcall & Soneira (1980). The open circles are the estimated 
number of galaxies on the NGC 6702 /-band frame (see text). The 
continuous line is the best fit to the faint /-band galaxy counts 
estimated by Woods & Fahlman (1997). At bright magnitudes 
the estimated number of galaxies detected on the NGC 6702 I- 
band frame are in fair agreement with those predicted by Woods 
& Fahlman (1997). At faint magnitudes (/ > 22.5 mag) incom- 
pleteness is affecting the galaxy counts from our study. 



length baseline of V — I compared to the B — I colour, pro- 
viding less metallicity leverage. The expected difference be- 
tween the peaks in V — I is smaller (i.e. A(V — I) ~ 0.3 mag; 
see next section) than that in the B — I histogram and there- 
fore random photometric errors can more easily smear out 
any bimodality. Contamination of the GC candidate sam- 
ple by stars and galaxies may also affect the significance of 
the detection of bimodality in the V — I colour distribution. 
This is demonstrated in Figure [l^, plotting the V — I colour 
distribution of contaminating stars and galaxies (after being 
selected on the basis of their magnitudes and colours). 



4.2 Age of the NGC 6702 Globular Clusters 

To constrain the age and metallicity of NGC 6702 GC sys- 
tem we employ a method similar to that developed by 
Whitmore et al. (1997) that uses the colour-magnitude rela- 
tion of GCs. The basic assumptions of the method are dis- 
cussed in detail by Whitmore et al. (1997). In brief, it is as- 
sumed that the GCs comprise two distinct sub-populations 
of different age and metallicity but with the same initial 
mass functions. An additional assumption is that both sub- 
populations have similar numbers of GCs. The older popula- 
tion is taken to have properties similar to those of the Milky 
Way GCs and therefore consists of old (15 Gyrs) metal-poor 
([Fe/H] = —1.5) GCs. The younger population is formed at 



Figure 12. B — I colour distribution for the candidate globular 
clusters, selected on the basis of their colours (see Figures g and ^) 
and magnitudes. Also shown is the colour distribution of Galactic 
GCs (hatched histograms). The blue peak of the NGC 6702 GCs 
lies close to the Milky Way GC distribution. 
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Figure 13. B — I histogram for the candidate globular clusters. 
The filled histogram is the distribution of contaminating stars pre- 
dicted from the models of Bahcall & Soneira (1980). The hatched 
histogram is the predicted colour distribution of the contaminat- 
ing galaxies within the candidate GC sample. 
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Figure 14. V — I colour distribution for the candidate globular 
clusters, selected on the basis of their colours and magnitudes. 
Also shown are the colour distributions of Galactic GCs (hatched 
histograms) . 




Figure 15. V — I histogram for the candidate globular clusters. 
The filled histogram is the distribution of contaminating stars pre- 
dicted from the models of Bahcall & Soneira (1980). The hatched 
histogram is the predicted colour distribution of the contaminat- 
ing galaxies within the candidate GC sample. 



a later stage of the evolution of the galaxy, through either a 
merger event (e.g. Ashman & Zepf 1992) or a secondary in 
situ star- formation burst (Forbes et al. 1997). The young 
GCs, assumed to be metal-rich, are initially bluer than 
the old metal-poor ones but become redder after about 1- 
2 Gyrs. Similarly, the newly formed GCs are much brighter 
than the pre-existing metal-poor population but become 
fainter at later stages as they evolve passively. Therefore, 
the difference between the colours and magnitudes of the 
two GC populations can be used to constrain the age of the 
new GCs. 

Here we employ the models of Worthey (1994) to predict 
the colour and magnitude difference of a young GC popula- 
tion relative to an old (15 Gyrs) metal-poor ([Fe/H] = —1.5) 
population. The predicted colours of this old stellar popu- 
lation are V - I = 0.921 and B - I = 1.564. The model 
predictions for the B — I colour against both B and / mag- 
nitudes are shown in Figures |l^ and |l^ respectively. For 
NGC6702 the difference A(B - I) between the colours of 
the two GC populations is taken as the difference between 
the peaks of the bimodal B — I distribution in Figure |l2[ 
The blue peak in that figure is assumed to correspond to the 
old metal-poor population. This is a reasonable assumption 
since the blue peak overlaps with the Milky Way GCs fairly 
well. The peak colour difference from the KMM statistical 
test is A(B — 7) = 0.64 ± 0.10, where the error is estimated 
assuming an uncertainty of 0.07 mag in defining the peaks 
of the bimodal distribution. 

The B and /-band magnitude difference, AB and AI, 
between the two GC populations is estimated following 
the method of Whitmore et al. (1997) by taking the off- 
set between the 10th brightest cluster of the blue and red 
subpopulations. The blue and red subpopulations are de- 
fined to have colours in the range 1 < B — I < 1.8 and 
1.8 < B — I < 3 respectively. The uncertainties in AB and 
AI axe estimating using the bootstrap resampling technique. 
We find AB = -0.05 ± 0.29 and AI = -0.84 ± 0.27, i.e. the 
young GCs are brighter than the old ones. 

The results are shown in Figures ^ and j^, and sug- 
gest that NGC 6702 GCs have solar to supersolar metallicity 
and an age in the range 2-5 Gyrs. Also we note that despite 
the errorbars, in both figures the data point lies close to the 
age of 3 Gyrs and [Fe/H]=0.5 model prediction. Moreover, 
although there are uncertainties in the calculation of the age 
of the newly formed GCs, our estimate is in fair agreement 
with the galaxy's spectroscopically estimated age and metal- 
licity of « 2 Gyrs and [Fe/H] w 0.5 respectively (Terlevich 
& Forbes 2000). 

To assess the sensitivity of the age and metallicity es- 
timates to the adopted single stellar population model of 
Worthey, we repeat our calculations using the models de- 
veloped by Bruzual & Chariot (1996) and Vazdekis et al. 
(1996) with a Salpeter IMF. Both models predict an age 
of about 2 Gyrs for the recently formed GC population. The 
predicted metallicity of [Fe/H] > 0.5 is higher than that from 
the Worthey model. Despite, the discrepancies, it is clear 
that all three models used here predict that the blue GCs 
sub-population is relatively young with a supersolar metal- 
licity. 

Additionally, assuming an age of 3 Gyrs and [Fe/H]=0.5 
for the red GCs, the Worthey (1994) models predict a V — 
I colour difference between the two GC subpopulations of 
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Figure 16. A(B — I) against AB diagram for the Worthey (1994) 
models. Each curve corresponds to a different mctallicity marked 
on the right of the curve. The crosses and the numbers on the 
curves indicate different ages in Gyrs. A(B — I) is defined as 
the colour difference between the red and the blue old (15 Gyrs) 
metal-poor ([Fe/H] = —1.5) GC populations. AB is defined as 
the magnitude difference of the 10th brightest globular cluster 
of the red and blue populations. The dot corresponds to the 
NGC6702 GCs (A(B - I) = 0.64 ± 0.10, AB = -0.05 ± 0.29), 
while the rectangle shows the la uncertainty region. 
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Figure 17. A(B — I) against AI diagram for the Worthey (1994) 
models. Each curve corresponds to a different mctallicity marked 
on the right of the curve. The crosses and the numbers on the 
curves indicate different ages in Gyrs. A(B — I) is defined as 
the colour difference between the red and the blue old (15 Gyrs) 
metal-poor ([Fe/H] = —1.5) GC populations. AI is denned as the 
magnitude difference of the 10th brightest globular cluster of the 
red and blue populations. The dot corresponds to the NGC 6702 
GCs (A(B - I) = 0.64 ± 0.10, AI = -0.84 ± 0.27) while the 
rectangle shows the la uncertainty region. 



A(V — I) ~ 0.3 mag. As discussed in the previous section 
such a small colour difference is difficult to detect with the 
present dataset. 

4.3 The surface density and specific frequency 

We study the surface density profile, p(r), of the NGC 6702 
GCs using all the sources detected in the B-band to B = 
26.1 mag (total of 1237), the 80% completeness limit for 
point sources. We use the B-band detections instead of the 
colour selected GC sample because the selection criteria (i.e. 
completeness limits) are easier to quantify and do not de- 
pend on the GC colours. The sample is binned into elliptical 
rings centred on NGC 6702 and having the same mean ellip- 
ticity and position angle as the galaxy (e = 0.2, PA=145°). 
The counts in each bin are corrected for geometric incom- 
pleteness due to foreground saturated stars, bad columns 
and the limited size of the field-of-view. Additionally, we 
mask out the central saturated part of NGC 6702 for radii 
r eq < 21". The correction for magnitude incompleteness to 
the counts in individual bins was found to be small (~ 3%) 
and was ignored. The density profile is plotted in Figure 
|l^(a). The background surface density level, ptg, was deter- 
mined by combining the counts within the outermost three 
bins lying at r eq > 120" (corresponding to galactocentric 
distances > 37 kpc) where the density profile reaches a con- 
stant level. We find p bg = (1.28 ± 0.05) x 10 5 sources deg" 2 



(assuming Poisson statistics) or 1114 ± 44 sources over the 
LRIS field-of-view of 31.2arcmin 2 (after applying geomet- 
rical corrections). The total number of objects to B = 
26.1 mag is 1237 and therefore the number of candidate 
GC is estimated to be 123±56 (Poisson statistics). In sec- 
tion [I we found a total of 151 candidate GCs, selected 
on the basis of their colours and magnitudes in the range 
22.5 < B < 26.1 mag. This is in fair agreement with the 
number of candidate GCs estimated here. 

Figure |l^(b) shows the candidate GC density profile, 
estimated by statistically subtracting the constant back- 
ground surface density level from the density profile of all 
the sources detected on the B-band frame to B = 26.1 mag. 
Assuming a power law density profile of the form pocr™ for 
the candidate GCs, we find a = —1.15 ± 0.52. This slope 
for the GC system of NGC 6702 is consistent with that for 
other galaxies of a similar magnitude (Forbes et al. 1997). 
Also shown in Figure [l8](b) is the B-band surface bright- 
ness profile of NGC 6702. It is clear that the starlight of 
NGC 6702 drops more steeply than the GC density profile 
as is often seen in other galaxies (Harris 1991; Durrell et 
al. 1996). This is often attributed to further dissipative col- 
lapse of the host galaxy after the older GCs had formed and 
before the halo stars formed (e.g. Harris 1986). The radial 
distributions of the red and blue subpopulations of the GC 
candidates selected on the basis of their colours and mag- 
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Figure 18. (a) Left panel: the surface density profile of all the sources detected in the NGC 6702 _B-band frame. The continuous lines 
is the mean background surface density, estimated by combining the counts in the bins at galactocentric distances > 120". The dashed 
lines are the errors around the mean assuming Poisson statistics, (b) Right panel: the surface density profile of the globular cluster 
candidates (open circles), derived by statistically subtracting the background surface density level from the density profile of the left 
panel. The best fit power law to the observed profile (continuous line) has the form p <x r — 1 - 1 s±0.52^ rp^g filled circles are the B-band 
surface brightness profile of NGC 6702 in arbitrary units. The starlight is more concentrated than the GC system. 



nitudes have also been investigated and have been found to 
be consistent within the uncertainties. 

The globular cluster specific frequency Sn is defined 
(Harris 1991) 



S N =N T xW 0A * (Mv+15 \ 



(1) 



where Nt is the total number of globular clusters in the 
galaxy and Mv is the host galaxy absolute l/-band lumi- 
nosity. To calculate Nt we assume that the B-band lumi- 
nosity function of NGC 6702 GCs is a Gaussian with the 
same parameters as that of the Milky Way GCs {Ms, V eak = 
—6.50 mag and o\b = 1.13). We note that differences in the 
peak B-band magnitude of the blue and red G C su bpopula- 
tion luminosity function are small (see section 4.2). At the 



distance of NGC 6702 MB, pe ak corresponds to B=27.5mag. 
Previously, we found that the number of GCs to the lim- 
iting magnitude B=26 mag is 123±56. Therefore integrat- 
ing over the entire area of the luminosity function we find 
Nt = 1776 ± 926. This implies Sn = 2.3 ± 1.1, assuming 
M v = -21.9 mag for NGC 6702. 

It is an interesting exercise to predict the evolution of 
Sn with time to investigate whether NGC 6702 will resem- 



ble normal ellipticals within a few Gyrs. For this calculation 
we assume that the total number of GCs is conserved and 
that the change in Sn is due to passive evolution of the 
host galaxy. Furthermore, we assume that during the most 
recent major star-formation event, that took place about 
2 Gyrs ago, 10% of the NGC 6702 mass (see section || below) 
has been converted into metal-rich ([Fe/H]=+0.5) stars. Af- 
ter 10 Gyrs this new stellar population will fade in the V- 
band by AV = 0.18 mag (Worthey 1994), corresponding to 
Sn = 2.7. Harris, Harris & McLaughlin (1998) found a mean 
of Sn = 5.4 ± 0.3 for large ellipticals in the Virgo, Fornax 
and the NGC 5846 group (excluding cD galaxies) . However, 
exclusion of Bright Cluster Galaxies (BCGs) results in a 
mean Sn = 3.3 ± 0.2, in better agreement, although some- 
what larger, with the estimated Sn ~ 2.7 for NGC 6702 
after 10 Gyrs. However, the estimated Sn evolution is sensi- 
tive to the galaxy mass fraction involved in the most recent 
major star-formation event, which in turn is uncertain and 
difficult to constraint. We calculate an upper limit for the 
NGC 6702 Sn after 10 Gyrs by assuming the extreme sce- 
nario where as much as 50% of the galaxy mass has been 
turned into stars. We estimate a fading of AV = 0.82 mag 
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NGC 6702 


NGC 1700 


Type 


E3 


E4 


iVJ -V \ L11CX, S) 


-21.9 


-21.5 


Environment 


Field 


Group 


(V/a)* 1 


0.21 


0.38 


Spectroscopic 






age 2 (Gyr) 


2 


2.3 


[Fe/H] 2 


0.5 


>0.5 


[Mg/Fe] 2 


0.12 


0.13 


Fundamental Plane 






residual 3 


-0.20 


-0.37 


o4 
°N 


2.3 ± 1.1 


1 4+ 1 - 



1 Bender et al. 1992 

2 Terlevich & Forbes 2000 

3 Forbcs, Ponman & Brown 1998 

4 The S N for NGC 1700 from Brown et al. 2000 



Table 2. Comparison between the properties of NGC 6702 and 
NGC 1700 (H = 75kms- 1 Mpc" 1 ). 



(Worthey 1994) within 10 Gyrs, corresponding to an upper 
limit of S N = 4.9. Therefore, we conclude that the NGC 6702 
Sn, although lower than that of bright cluster ellipticals, it 
is likely to be similar to that of 'normal' (old) present-day 
ellipticals (excluding BCGs) within the next 10 Gyrs. The 
calculations above are also consistent with suggestions that 
galaxy interactions and mergers although unable to account 
for large Sn galaxies, may explain the low values of Sn ~ 2 
in relatively isolated ellipticals (Elmegreen 1999). 



5 DISCUSSION 

A comparison between the properties of NGC 6702 and the 
group proto-elliptical NGC 1700 (Brown et al. 2000; Whit- 
more et al. 1997), is shown in Table 2. It is clear that there 
are quite a few similarities between the two galaxies. Both of 
them have similar optical luminosities and deviate from the 
Fundamental Plane in the sense of having younger ages on 
average (Forbes, Ponman & Brown 1998) . This is confirmed 
by their central spectroscopic ages of ~2 Gyrs (Terlevich & 
Forbes 2000). The central stars appear to be of high metal- 
licity ([Fe/H] ~0.5) with supersolar [Mg/Fe] abundance. Su- 
persolar abundances have been the source of much debate 
in the literature but probably indicate either a rapid star 
formation time-scale or an IMF skewed towards massive 
stars. Both galaxies have radial velocity to velocity disper- 
sion ratios, (V/a)* , indicating that they are anisotropic ro- 
tators. Bender et al. (1992) have suggested that the last 
major merger for anisotropic galaxies involved relatively lit- 
tle gas, i.e. the progenitors were largely stellar. In the case 
of NGC 6702 and NGC 1700 there must have been sufficient 
gas to support the central starburst and associated cluster 
formation. The similarities above also suggest a similar for- 
mation mechanism for these galaxies. 

Moreover, our analysis above indicates the presence of 
two GC populations in NGC 6702 with different metallici- 
ties and ages. It is an interesting exercise to estimate the 
formation efficiency of the young (red; 1.8 < B — I < 3) 
GC subpopulation, assuming they formed during the most 



recent star-formation burst, and compare it to that of stars 
formed in the same burst. 

A robust estimator of the GC formation efficiency is 
the ratio of the total GC mass to the total baryonic galaxy 
mass. For this calculation we assume a typical GC mass of 
3 x 10 5 M Q (Harris, Harris & McLaughlin 1998) and a mass- 
to-light ratio for NGC 6702 (M/L) v = 8M /L e , typical 
to that of ellipticals (Faber & Gallagher 1979; Binney & 
Tremaine 1987). We have identified a total of 88 red GCs 
(1.8 < B - I < 3.0) to the limit B » 26.0 mag, assumed 
to be the newly formed metal-rich GCs of NGC 6702. We 
adopt a Gaussian luminosity function for this GC subpop- 
ulation with parameters similar to those used in section 4.3 
(MB, P eak ~ —6.50 mag and o\b = 1.13). This is a reasonable 
assumption, since the blue and red GC subpopulations of 
NGC 6702 in section 4.2 have a similar B-band magnitude 
distribution within few tenths of a magnitude. Integrating 
over the entire luminosity function, assuming a complete- 
ness limit for our sample B ~ 26.0 mag, we estimate a total 
of 970 red GCs. Thus, the total mass of the red GCs is 
Mac = 970 x 3 x 10 s M @ = 2.9 x 10 8 M Q . The total bary- 
onic mass in stars of NGC 6702, estimated using its l/-band 
luminosity and the assumed mass-to-light ratio, is found to 
be M G = 3.6 x 10 11 M Q . Therefore, we estimate that 0.08% 
of the NGC 6702 mass has been turned into GCs during the 
most recent star-formation event, which is lower compared 
to recent estimates for the GC formation efficiency of 0.25% 
(Harris, Harris & McLaughlin 1998; McLaughlin 1999). It 
should be noted that this calculation assumes that the es- 
timated baryonic mass in stars represents the total mass of 
the system, i.e. most of initial gas has been used for star- 
formation. 

Elmegreen et al. (1999), argue that GCs represent about 
1% of the total galaxy star- formation. Therefore, for the 
GC formation efficiency of 0.08% we expect that « 8% of 
the total NGC 6702 mass has turned into field stars (not in 
bound clusters) during the most recent star- formation event. 
Is our estimate of the mass in new stars (i.e. 8% of the to- 
tal) consistent with the spectral line indices in the central 
r e /8 region (Gonzalez et al. 1993)? If we assume that the 
new stars formed with a metallicity of [Fe/H]=0.5, and that 
the old stars (representing 92% of the central stellar mass) 
have solar metallicity and an age of 15 Gyrs, then the new 
stars would have to be slightly younger than 1 Gyr old. As- 
suming a time since the starburst of « 2 Gyrs and applying 
that to the new stars, their mass fraction would have to be 
increased to « 15%. Furthermore, if the galaxy reveals ev- 
idence for a young stellar population outside of the central 
region this will also point towards a large mass fraction (see 
also Trager 1997). Clearly there are large uncertainties in- 
volved in the above arguments, but there is some agreement 
between the mass fraction of young stars estimated inde- 
pendently using spectroscopic line indices and the proper- 
ties of the NGC 6702 GC system. Moreover, it seems clear 
that a large fraction of new stars formed in NGC 6702 a 
few Gyrs ago. In the merger context, it implies a high gas 
fraction for the progenitors and hence a highly dissipational 
event. This is contrary to the claims of Bender et al. (1992) 
that anisotropic galaxies (such as NGC 6702 and NGC 1700) 
formed by largely dissipationless mergers. 

Next, we compare the fraction of galaxy starlight in 
newly formed dense clusters relative to that of young field 
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stars. To do that we add the luminosities of the red GCs 
(total of 88; Lgc,b = 4.9 x 10 7 Lq) and then divide by 
the luminosity of the young stellar population of NGC 6702. 
We note that estimating the total GC luminosity by extrap- 
olating to the GC luminosity function will not modify our 
results, since the main contribution to the luminosity is from 
bright systems. Therefore, adding the red GC luminosities is 
a simple approximation that is sufficient for the purposes of 
the present study. Additionally, in our calculation we assume 
that 10% by mass of the NGC 6702 stars are young (2 Gyrs) 
and metal rich ([Fe/H]=0.5), while the remaining 90% of 
the mass is in old (15 Gyrs) metal poor ([Fe/H]=-1.5) stars. 
Using the mass-to-light ratios of the above stellar popula- 
tions predicted by the models of Worthey (1994) we estimate 
that 15% of the total NGC 6702 B-band light originates in 
newly formed stars, i.e. I/*,s = 2.9 x 10 9 Lq. Consequently, 
we estimate Lgc,b /L+,b ~ 0.02 , i.e. 2% of the young star 
B-band luminosity originates in bound clusters formed dur- 
ing the most recent star-formation event. This is much lower 
than the estimates in recent mergers which have shown that 
as much as 20% of the blue galaxy light originates from 
GCs (Zepf et al. 1999; Meurer 1995). However, the main 
contribution to the integrated starlight in bound clusters in 
these systems is from young luminous clusters, the nature 
of which is under debate. In particular, it is still unclear 
whether these objects are proto-globular clusters that will 
evolve into bona-fide GCs (Brodie et al. 1998) or are lu- 
minous open clusters that will be destroyed by dynamical 
evolution (van den Bergh 1995). Indeed, in older merger re- 
mants, like NGC 7252, the fraction of light in bound clusters 
is somewhat lower than 20% (Zepf et al. 1999). Moreover, 
van den Bergh (1994) noted that GCs account for about 2% 
of the total halo luminosity of both the Milky Way and the 
LMC, in good agreement with our estimate for NGC 6702. 

Clearly, the young elliptical galaxy NGC 6702 (and also 
NGC 1700) is at an interesting evolutionary stage for further 
study of its GC system. Direct measurement of individual 
GC metallicities with 8- 10m class telescopes would confirm 
(or refute) the inferred supersolar metallicities for the new 
GCs and hence their age estimates. We n ote t hat the mag- 
nitude differences AB and Al in section 4.2, used to con- 



strain the GC age and metallicity, are hard to define and 
would benefit from deeper imaging. Deeper imaging stud- 
ies would also help determine the GC luminosity function 
(e.g. Gaussian or power law) to explore possible destruction 
mechanisms and would allow a detailed study of the radial 
distribution of the GC subpopulations. 



6 CONCLUSIONS 

We have carried out B, V and /-band photometric obser- 
vations of the elliptical galaxy NGC 6702 to explore its GC 
system. The GC candidates are selected on the basis of their 
colour, magnitude and morphology. We find strong evidence 
that the B — I colour distribution of GC candidates is bi- 
modal with the blue peak having a colour similar to that of 
the old metal-poor Galactic GCs. Assuming that the blue 
GC subpopulation is indeed old and metal-poor, we esti- 
mate that the red GCs have an age of 2-5 Gyrs and super- 
solar metallicity. Despite the large uncertainties, this is in 
reasonable agreement with the spectroscopically estimated 



age and metallcity of about 2 Gyrs and 0.5 respectively for 
the galaxy. Thus the red GCs in NGC 6702 are consistent 
with having formed during the most recent starburst, likely 
to have been triggered by a merger or an accretion event, a 
few Gyrs earlier. Additionally, we find a specific frequency 
for the NGC 6702 GC system of S N = 2.3 ± 1.1. We pre- 
dict that passive evolution of the NGC 6702 stellar popu- 
lation formed during the most recent star-formation event 
will increase the NGC 6702 specific frequency to « 2.7 within 
10 Gyrs, in reasonable agreement with the mean for present- 
day ellipticals, excluding the largest cluster ellipticals. Fi- 
nally, using the NGC 6702 GC properties and making rea- 
sonable assumptions about the GC formation efficiency, we 
estimate the mass fraction of newly formed stars during the 
most recent star-formation event. We find reasonable agree- 
ment with the mass fraction estimated using stellar popu- 
lation synthesis models to interpret the line indices of the 
underlying galaxy starlight. 
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